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Abstract 

In the presence of weak time harmonic electromagnetic fields, Superconducting Quantum Interference 
Filters (SQIFs) show the typical behavior of non linear mixers. The SQIFs are manufactured from high- 
T c grain boundary Josephson junctions and operated in active microcooler. The dependence of dc voltage 
output Vdc vs. static external magnetic field B is non-periodic and consists of a well pronounced unique dip 
at zero field, with marginal side modulations at higher fields. We have successfully exploited the parabolic 
shape of the voltage dip around B = to mix quadratically two external time harmonic rf-signals, at 
frequencies /i and f2 below the Josephson frequency /j, and detect the corresponding mixing signal at 
I fi — /2 I- When the mixing takes place on the SQIF current- voltage characteristics the component at 
2/2 — fx is present. The experiments suggest potential applications of a SQIF as a non-linear mixing 
device, capable to operate at frequencies from dc to few GHz with a large dynamic range. 
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1. INTRODUCTION 



Superconducting Quantum Interference Filters (SQIFs) are serial (or parallel) arrays of Joseph- 
son junctions with an irregular grating structure, i.e. a specially selected distribution of the loop 
areas, such that the phase sensitive superposition of the individual loop modulations leads to a 
dc voltage response V dc vs. static magnetic field B [V dc (B)- curve] that has a unique dip around 
B — 0, with a vanishing modulation everywhere else [1, 2]. The detailed shape of the voltage 
output Vd c (-B), and in particular the dip width AB and its voltage span AV, can be engineered by 
choosing both suitable areas of the loop sequence and electric parameters of the Josephson junc- 
tions. However, the SQIF performance is not degraded by the large spread of parameters typically 
present in high T c -junctions. The sensitivity of a SQIF in response to an applied magnetic field 
is determined by the voltage-to-magnetic field transfer factor Vb = max(dVd c /dB). Employing 
various flux focusing structures together with a SQIF it is possible to significantly enhance the 
transfer factor [3, 4]. A flux focusing structure integrated together with a SQIF array is in the 
following referred to as a SQIF-sensor [for a schematic view see Fig. 1(c)] . 

In the resistive mode, the SQIF dc voltage output V dc (Ib, B) is set by the bias current 1^ and 
by the static field B, and it is the time average of a fast signal V(It,,B,t) oc d<p/dt, where 
if is the Josephson phase across the JJ. The main harmonic of V(Ib, B,t) has frequency fj = 
Vdc(^b, B) I $ (the second Josephson relation) [7]. 

Let us consider the case when in addition to the static magnetic field B applied perpendicular to 
the SQIF loops, there is a weak radiofrequency magnetic field b r f (t), such that B(t) = B + b r { (t). 
The time harmonic magnetic field is accompanied by a time harmonic electric field, E{t) = E + 
e r f (t). If we assume the electric field to be oriented parallel to the bias leads, the corresponding 
fluctuations induced in the bias current can be written as I(t) = lb + iri{t), which are naturally 
dispersive due to the frequency dependence of the surface impedance in a superconductor. Thus, 
if the fields vary slowly in comparison with fj, the SQIF voltage can be written as 



where <9i b [<9B] denotes the derivative with respect to 1^ [B], and () the time average over the 
time scale set by fj 1 , while the slow time dependence remains. In Eq. (1) we have denoted as 




bias coupling 



flux coupling 



V dc (I h ,B) + Mt)d lh + b T{ (t)d B )V dc (I h ,B) + -Mt)d lh + b T f(t)d B yV dc {I h ,B) + ... 
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bias coupling the dependence of the SQIF voltage on the fluctuating induced rf-current, and as 
flux coupling the dependence on the fluctuating magnetic flux. For 6 r f (t) = 6 r f cos(2ir ft), with 
amplitude b l{ <C AB and frequency / < 4^, Eq. (1) describes a superposition of the dc voltage 
output Vdc(/b, B) with slowly varying terms in / (voltage outputs at frequency /, due to bias 
and flux variations, with pre-factors proportional respectively to the amplitudes z r f and 6 r f); plus 
rf terms in 2f, etc... Thus, the amplitude of the first harmonic of the Fourier spectrum of the 
output signal contains information about the slopes #BVd c |i b and <9i b Vd c |B> while the amplitude of 
the second harmonic is proportional to the curvature in the 9gVd c |i b and in the d? Vdcb • 

The Fourier transform of V(Ib + i r f (t) , B + b T f(t) , t) is the spectral voltage output V{Ik, B, /). 

~ 2 

Using a spectrum analyzer, the power spectrum V ( Jb, B, f) can be detected at frequencies 
around the center frequency /, as a function of the sweeping parameter which can be either 1^ (at 
a constant value of the static field B), or as a function of B (at constant value of the bias current 
Jb)- If the amplitude of the incident rf signal is sufficiently small, b ri <C AB, the spectral voltage 
output will be proportional to $BVdc|i b and to <9[ b Vd c |B: it will be maximum in field regions of 
maximum slope, while in regions of a flat Vd c (Ib, B) curves, it will tend to zero. 

Quadratic mixing is expected for two tone experiments, where the incident rf signal is a su- 
perposition of two time harmonic signals: 6 rf (t) = 6 rfl cos (2nfi t) + b ri2 cos (2irf 2 1) (we assume 
/i < $2 < 2 )• Taking into account the parabolic shape of the Vd c (B)- curve around B = 0, one 
expects in that region a larger amplitude for the quadratic mixing signals at frequencies f 2 — fi, 
2/i, f2 + fi and 2f 2 , respectively, compared to the field regions of zero curvature of Vdc(-B)- 
characteristics. In a similar manner, in the non linear regions of the Vdc(-^b)- curve, quadratic 
mixing can also take place. However, in the case of the Vdc(-^b) -curve the spectral voltage output 
contains higher order mixing products, due to higher order non linearities present in the character- 
istics. 

Recently, we have demonstrated quadratic mixing of radio frequency signals using SQIF- 
sensors [6]: the external rf signals couple through the time dependent magnetic flux threading 
the loops of the interferometer, giving rise to a beat-pattern of the original high frequency voltage 
output (at the Josephson frequency), set by I h and B. The quadratic mixing takes place due to exact 
parabolic shape of the V<± C (B)- curve, and the mixed signal is detected at frequencies | fx ± / 2 |. 
Due to limited bandwidth of the set-up, the incident signal frequencies were always chosen such 
that the quadratically mixed components were at a few hundred MHz. The experiments were per- 
formed in active micro-coolers, without significant degradation of the SQIF-sensors, even when 
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the magnetic shield was removed. 

In this paper, we present experiments in which the mixing of radiofrequency signals on the 
voltage output of the SQIF is extended to larger frequencies. The mixing is studied in two ways: 
either using as sweeping parameter the static magnetic field B {i.e. mixing made on the V& C (B)- 
curve, at a constant value of the bias current 1^), or the bias current I h (i.e. mixing made on the 
Vdc(-fb)- curve, at B = 0). As expected, also when the incident fields have larger frequencies the 
detected spectral voltage output of our SQIF sensor shows quadratic components, with maximum 
amplitude in the regions of maximum curvature of the Vd c (B) curve. In the regions of non linear 
Vdc(-fb)- curve, the quadratic mixing effect takes place as well. This can be explained by the 
dependence of the SQIF output voltage on the time varying electric field [through the bias current, 
as described by Eq.(l)]. Only when mixing on the Vd c (Ib)- curve, high order mixing components 
are also observed. These are never present when mixing on the Vdc(-B)- curve - as long as small 
amplitude signals are used (6 rf <C AB). 

2. THE SETUP 

Experiments have been performed with serial SQIFs of 21 1 loops, consisting of YBa 2 Cu307_ x 
grain boundary Josephson junctions grown on 24°- oriented bicrystal MgO substrates[ ]. The 
Josephson junctions are designed with a width of 2 /im, the YBa 2 Cu 3 07_ x layer being 130 nm 
thick, so that the resulting junction critical current density is J c « 23kA/cm 2 , at T = 77 K. The 
serial loops have an area distribution ranging between 38/im 2 and 210 /im 2 . In order to enhance the 
coupling of the SQIF loops to the magnetic flux, we fabricated a SQIF-sensor with a flux focusing 
structure, i.e. a superconducting ring placed on one side of the grain boundary. Figure 1(a) is an 
optical microscope image of the sample. The focusing loop is effectively a split loop design[13], 
consisting of 10 equidistant parallel thin loops, one aligned inside the other; a fragment of the split 
loop along the SQIF is enlarged in Fig. 1(b). Figure 1(c) shows a sketch of the SQIF-sensor with 
the serial loops of the interferometer along the grain boundary, and the inductively coupled flux 
focusing YBa 2 Cu30 7 „ ;r ring. When the pick-up ring is split in various parallel loops, we have 
obtained an enhanced sensitivity of the SQIF to the magnetic field (larger transfer factor) with 
respect to the case of a single broad ring. Calculations show [13] that the split ring design leads, 
in average, to a larger current density compared to a single ring of the same cross section. 

The static magnetic field B is applied via a multi-turn coil placed inside the cooler. The rf 
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(c) B + b rf 

FIG. 1: (a) Optical microscope image of the chip with the SQIF placed along the grain boundary, and the 
superconducting ring inductively coupled to the SQIF; (b) zoom on the ring made of 10 equidistant parallel 
loops, one aligned inside the other. All superconducting parts, except the regions across the junctions, are 
covered by gold, (c) Sketch, not in scale, of the SQIF-sensor, indicating the grain boundary (dashed line), 
the direction of the bias current and of the magnetic field. 

magnetic field b T { (i) is applied via a 50 fi-loop antenna (primary antenna). Experiments are made 
either with a mu-metal shield - in this case the rf antenna was placed inside the shield at a distance 
of about 5 cm from the chip -, or with an unshielded cryocooler in open space. 

Samples have been operated in a Stirling microcooler [8], at temperatures from 55 to 82 K. 
Active cooling does not degrade the SQIF performance, and offers numerous advantages, such 
as the possibility to set stable temperatures over a large temperature range, quick thermal cycles, 
compact and portable setups [9, 10]. To provide transparency to the incoming rf fields, a glass 
chamber surrounds the microcooler cold head. The sample is anchored to the cold head only from 
one side, in order to keep the focusing ring as much as possible free of metallic ground plane. On 
the cold head, an rf amplifier is also mounted, designed for cryogenic applications with a gain 
of about 40 dB, for the bandwidth 0.04 — 6 GHz. At room temperature and at the typical bias 
values, we have measured a noise figure of the amplifier of about 1 dB, in a frequency range from 
100 MHz to 1 GHz. The amplifier input resistance is 5 kfi. The SQIF is coupled to the amplifier 
through a coplanar line, after which a bias-Tee is placed. The bias-Tee splits the ac and the dc 
parts of the SQIF signal: the ac part, through a capacitor, follows the input line to the first stage 
of the cold amplifier; the dc part (bias input & voltage output), through a conical inductor, goes 
to the room temperature electronics, via twisted wires. The amplifier output is carried outside the 
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cooler via an SMA feed-through and flexible coaxial lines, and detected by a spectrum analyzer. 

The rf measurements are performed in the following way: the spectrum analyzer is operated in 
Zero Span Mode, i.e. as a narrowband receiver tuned at the central frequency f (equal to the fre- 
quency at which the detection of the SQIF spectral voltage output is made), and with a bandwidth 
around f set by the Resolution Bandwidth (ResBW). To this mode, it corresponds an analog out- 
put proportional to the maximum amplitude of the signal at f . This output voltage is recorded by 
computer while slowly sweeping either the static magnetic field B (simultaneously, the V& C (B)- 
curve is acquired), or sweeping the static bias current J b (simultaneously, the V dc (I h )- curve is 
acquired). 

3. THE EXPERIMENTS 

3.1 Mixing on the SQIF V dc (B)- curve 

We first measure the Vd c (B) dependence of the SQIF at a constant 7b value, and check the 
dip symmetry with respect to zero magnetic field. There is an optimal temperature at which the 
voltage swing is a maximum, at a properly chosen bias current. At T « 73 K and at I h = 20 /xA, 
we measure a voltage span AV ~ 1160/iV and a transfer factor Vb ~ 6500 V/T. The SQIF 
normal resistance is R = 186 fl Successively, the rf is switched ON: a time harmonic signal 
which is the superposition of two rf signals with frequencies fi and f2 and equal amplitudes 
6 r fi and fe r f 2 is applied to the primary antenna, while the rf output of the SQIF is detected. A 
small amplitude of the incoming signals is required, so that the rf signal superimposed to the 
static field does not modulated the SQIF working point out of the dip and corrupt the effect; but 
it has to be high enough so that the SQIF spectral voltage output is above the noise level set by 
the resolution bandwidth (ResBW) of the spectrum analyzer. The analog output of the spectrum 
analyzer is recorded by computer while slowly sweeping the static field B (sweep frequency in 
the kHz range); simultaneously, the dc voltage output Vd c is acquired as a function of B. In the 
output spectra, we find components at f =\ fi — f 2 |, whose amplitude is maximal at the dip 
bottom and vanishes in the region of zero curvature of the Vd c (B) curve. Indeed, the amplitude of 
the second harmonics (the term in fi — f 2 ) is expected to be proportional to the curvature in the 
<9!^dc|i b - We have covered the frequency interval from few hundred MHz to few GHz for the two 
primary incident fields, keeping the difference frequency within the bandwidth of the rf amplifier, 
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for example at values f x — f 2 = 1.5GHz, but also 2 and 2.5 GHz. 

Figures 2 displays the typically measured curves. The incident field is a linear combination of 
two signals with frequencies /i = 3 GHz and f 2 = 1 GHz, and equal amplitudes. In the spectral 
voltage output of the SQIF, we find then a quadratic mixing signal at the difference frequency 
fo = 2 GHz, and in the Figure is displayed the spectral voltage output detected at the central 
frequency / , with ResBW equal to 30 kHz, as a function of the sweeping parameter B. The 
amplitude of the voltage output at f is maximal around B = 0, and vanishes as the dc working 
point approaches the dip slopes. 



140 - 



- 3.0 




FIG. 2: The two incident signals are set at f\ = 3.0 GHz and f 2 = 1.0 GHz. The rf voltage is measured 
with a ResBW of 30 kHz. 1^ = 20 /uA and T = 73K. Gray curve, relative to right axis: V dc vs. B ; black 
curve, rel. to left axis: rf voltage V (B, lb, /) detected at fo = f± — f 2 = 2 GHz. 



3.1 Mixing on the SQIF V dc (I h )- curve 

At B = 0, the V dc (I h ) is swept; the two incident rf signals are applied to the primary antenna, 
and the spectral voltage output of the SQIF sensor is detected. To the applied time harmonic 
magnetic field, there corresponds naturally a time harmonic electric field, which couples to the 
SQIF voltage output through the fluctuations induced in the bias current (if we assume the electric 
field to have a component parallel to the bias line). Thus, the spectral voltage output detected 
vs. Ib might also contain fingerprints of the incident signals, and eventually mix them. Figures 
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FIG. 3: The two incident signals are set at fi = 3.5 GHz and /2 = 2.2 GHz. The rf voltage is measured 
with a ResBW of 3 kHz. B = and T = 73.5 K. (a) Gray curve, relative to right axis: V^ c vs. 1^, ; 
black curve, rel. to left axis: rf voltage V (B, 1^, /) detected at / = f\ — ji = 1.3 GHz; (b) Gray 



curve, relative to right axis: V^c vs - lb ; black curve, rel. to left axis: rf voltage V (B, 1^, f) 
fo = 2/2-/1 = 0.9 GHz. 



detected at 



3(a-b) display the typically measured curves. The incident field is a linear combination of two 
signals with frequencies /i = 3.5 GHz and f 2 = 2.2 GHz, and equal amplitudes. In the spectral 
voltage output of the SQIF, we find the quadratic mixing component at the difference frequency 
/o = fi — h — 1-3 GHz. In Fig. 3(a) is reported the voltage output detected at f = 1.3 GHz, with 
a ResBW of 3 kHz, as a function of the sweeping parameter 1^. Also in this case, the amplitude 
of the voltage output at f\ — f 2 is found to be maximal around the region of maximal curvature 
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of the Vdc(-fb)- characteristics. In contrast to the case of the Vdc(-B)- curve, when sweeping the 
Vdc(-^b)- curve the spectral voltage output contains also a component at f = 2/ 2 — f\. Figure 3(b) 
shows the SQIF rf voltage output detected at f = 0.9 GHz. As expected, while the amplitude of 
the term in /i — f 2 is proportional to the curvature in the <9f Vd c |B, the amplitude of the mixing 
component at frequency 2/ 2 — fi is proportional to the absolute value of the 3-rd order derivative 
of the Vdc(Jb)- curve. 

We have also performed mixing experiments without magnetic shield, both on the Vd c (B)- and 
on the Vd c (Ib)- curve. The presence of environmental disturbances does not suppress the effect 
(although in this case the compensation field is different than zero), and the mixing components 
are clearly visible. 

CONCLUSIONS 

In conclusion, we have demonstrated quadratic mixing of radiofrequency signals using the 
voltage output of a SQIF, both on the Vd c (B)- and on the Vd c (/b)- curves. The SQIF sensor 
is capable to transfer an incident radiofrequency signal into a converted radiofrequency voltage 
signal: the external rf fields couple either through the time dependent magnetic flux threading the 
loops of the interferometer or through the bias current, giving rise to a beat-pattern of the original 
high frequency voltage output at the Josephson frequency. In the adiabatic approximation, the 
quadratic mixing of the incident rf signals takes place due to the exact parabolic shape of the dc 
voltage output [V dc (B)- curve] around B = 0. When mixing on the Vdc(-^b)- curve, in addition to 
the quadratic components we have found higher order mixing products. In all cases, the detected 
harmonics have an amplitude dependence proportional to the corresponding order derivative of the 
Vd c - curve with respect to the sweeping parameter. 

In the presented experiments, by means of a SQIF sensor we have obtained quadratic mixing of 
incident fields in a frequency range from few hundred MHz up to several GHz, without significant 
degradation of performance. This feature might open new possibilities for applications of SQIF 
sensors as broad band mixers. 

The authors are very thankful to R. IJsselsteijn to provide us with samples, and V. Schultze for 
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